Entropic effects in the self-assembly of open lattices from patchy particles 
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Open lattices are characterized by low volume-fraction arrangements of building blocks, low coor- 
dination number, and open spaces between building blocks. The self-assembly of these lattices faces 
the challenge of mechanical instability due to their open structures. We theoretically investigate the 
stabilizing effects of entropy in the self-assembly of open lattices from patchy particles. A prelimi- 
nary account of these findings and their comparison to experiment was presented recently [X. Mao, 
Q. Chen, S. Granick, Nat. Mater., 12, 217 (2013)]. We found that rotational entropy of patchy 
particles can provide mechanical stability to open lattices, whereas vibrational entropy of patchy 
particles can lower the free energy of open lattices and thus enables the selection of open lattices 
verses close-packed lattices which have the same potential energy. These effects open the door to 
significant simplifications of possible future designs of patchy-particles for open-lattice self-assembly. 

PACS numbers: 81.16.Dn, 65.40.gd, 82.70.Dd, 46.32. +x 
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I. INTRODUCTION 

Open structures, structures characterized by low 
volume-fraction arrangements of building blocks with 
open spaces between them, are ubiquitous in both natu- 
ral PHD] and engineered |8l-KL2| systems. Open-structure 
materials occur in a rich variety of morphologies^ and 
exhibit fascinating aspects in their mechanical pi 
transport [H,[l7]], and thermodynamical properties 
leading to important applications in engineering. For 
example, zeolite a natural aluminosilicate mineral, 
consisting open-lattice structure with pores of sizes at 
the nanometer scale, has been widely used as catalyst, 
molecular sieves, solar thermal collectors, etc., and has 
been shown to exhibit a number of unusual features in- 
cluding negative thermal expansion Q , negative Poissons 
ratio pj| Pil [2(| , and flexibility window pj| . 

Structural openness of a material often leads to struc- 
tural flexibility, which is a double-edged sword. On the 
one hand, open structures are often close to mechanical 
instability point and this poses great challenges to the 
fabrication of open-structure materials. On the other 
hand, the closeness to mechanical instability point also 
opens the door to interesting designs of novel materials 
with interesting applications such as precisely controlled, 
fast, and reversible deformations, etc. 

The mechanical instability of open-structure materials 
can be understood through a simple counting argument 
due to J. C. Maxwell in 1864 [2l|. The basic idea of this 
counting argument is that to determine whether a system 
is mechanically stable one can simply compare the total 
number of degrees of freedom and total number of 
constraints N c . If N d > N c + d(d+l)/2, where d(d+l)/2 
represent the degrees of freedom of the global translations 
and rotations of the whole system, this system does not 
have sufficient number of constraints to fix all the rela- 
tive- degrees of freedom and there exist deformations, i.e., 
ways to rearrange the different components of the system, 
without violating any of the constraints. These deforma- 
tions are called "floppy modes" and we use Nq to denote 



the number of these floppy modes. The point at which 
Nd and N c + d{d + l)/2 (assuming all constraints non- 
redundant) becomes equal defines a mechanical critical 
point, at which the system is at the onset of mechanical 
stability. This point is usually called the "isostatic" point. 
Above this point, the system exhibit mechanical stability 
and there are no floppy modes. For an actual system, 
constraints often take the form of elastic connections, 
e.g., contacts between colloidal particles, sections of semi- 
flexible polymers between crosslinks, struts in frames, or 
chemical bonds between atoms. Correspondingly, floppy 
modes take the form of zero-energy deformations, which 
lead to mechanical instabilities in some form. Other non- 
floppy modes in the system are allowed but will cost elas- 
tic energy. 

In particular, for periodic lattice structures consisting 
of point-particles with only nearest neighbor central-force 
interactions, this criterion of no floppy modes can be sim- 
plified by ignoring the global translational/rotational de- 
grees of freedom d(d+ l)/2 given that it is much smaller 
than Nd and N c in a large lattice. Thus one have the sim- 
ple equation for isostaticity in periodic lattices, z = 2d, 
where z is the coordination number, i.e., the number of 
connecting neighbors of each particle. At this point, each 
lattice site has just enough constraints to fix all its de- 
grees of freedom. In practice, lattices are of finite size, 
and particles on the boundary have fewer constraints 
than the ones in the bulk. As a result there exist a sub-ex- 
tensive number of floppy modes in a finite lattice in which 
z = 2d in the bulk. These floppy modes can either take 
the form of extended or surface floppy modes, depending 
on the geometry of the lattice [l3j . 

Periodic open-structure lattices, such as honeycomb, 
kagome lattices in two dimensions and diamond, py- 
rochlore, perovskite, diamond lattices in three dimen- 
sions, are typically either at or below the isostatic point, 
thus at finite size they exhibit floppy modes which lead 
to mechanical instability and makes them challenging to 
produce in practice, as we mentioned above. In existing 
open-structure materials, additional interactions beyond 
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nearest-neighbor central-force interactions provide more 
constraints which stabilize the open lattices. For exam- 
ple, in open-lattice crystals, covalent bonds with prefer- 
ence on bond-angle directions stabilizes spatially open 
structures such as the zeolite or crystobalite, whereas 
in bio-polymer networks the molecular structure of the 
polymers provides stiffness against bending and thus pro- 
vides mechanical stability. 

At the colloidal scale, these stabilization mechanisms 
are difficult to realize, because interactions are typically 
nearest-neighbor without bond-angle selection. Provided 
the wide predicted applications of open structures at this 
scale in fields such as catalysis and photonics p3 - |26| . 
many designs have been proposed to obtain open-lattices 
from self-assembly in the colloidal scale via additional 
potential energy terms such as further range interaction, 
three- or more-body interaction, DNA-functionalization, 
etc. p7l - [33l ]. Most of these designs are still difficult to 
realize given current experimental techniques on colloidal 
particles. 

In this Paper, we discuss a new mechanism of sta- 
bilizing open colloidal lattices through cntropic effects 
in patchy particles, and present in detail a systematic 
theory that captures this effect and provides guidelines 
to remarkably simple future designs of open-lattice self- 
assembly. We apply this theory to the model system of a 
two-dimensional self-assembly of triblock Janus particles 
which arc particles with attractive patches on the north- 
and south- polar caps [26|, HH , and arrive at free energies 
of different lattice geometries and the equilibrium phase 
diagram under pressure. This theory is readily generaliz- 
able to other systems of self-assembly of patchy particles. 

Using the Maxwell's counting argument discussed 
above, the two-dimensional kagome lattice is at the iso- 
static point. At finite size with free boundary condition, 
there exist a sub-extensive number of floppy modes. For 
the kagome lattice these floppy modes correspond to the 
rotations of triangles along straight lines in the lattice as 
shown in Fig.[TJ These modes are of zero potential-energy 
cost if there are only nearest-neighbor central-force in- 
teractions in the lattice, because to leading order they 
keep the length of all bonds unchanged. Finite ampli- 
tude floppy modes also exist in the kagome lattice and 
are discussed in Ref. [l3j . 

The tri-block patchy particles do not directly provide 
additional potential energy penalties to eliminate these 
floppy modes [H, ® • Patchy particles are characterized 
by chemical coatings ("patches") that cover finite frac- 
tions of their surfaces, and the coated parts of the surface 
exhibit interactions of very different nature from the un- 
covered parts. The interactions between these particles 
are typically very short-ranged, and thus due to the ex- 
tended sizes of the patches, these patchy particles often 
do not have a direct selection of bond-angle directions, as 
we will discuss in detail in Sec.|TTJ The potential energy of 
pair-wise interactions between these patchy particles are 
of step-function nature as relative particle-orientations 
are changed, rather than gradual dependence with a po- 




FIG. 1. The kagome lattice and one example of its floppy 
modes. The floppy-mode deformed kagome lattice is shown 
in solid lines, and the undeformed lattice is in dashed lines. 
The deformations only occur in a horizontal line in which the 
pointing-up and pointing-down triangles are rotated. In other 
parts of the lattice the deformed and the undeformed lattices 
are on top of each other. 



tential energy minimum at certain angle. 

For the case of tri-block Janus particles discussed in 
Rcfs. 26, 34], two particles experience a short-ranged mu- 
tual attraction when they are facing each other in their 
coated attractive patch, and strongly-screened electro- 
static repulsion (which can be modeled as hard-sphere 
repulsion) otherwise. The pair- wise interaction potential 
changes abruptly at the boundary of the patches, and is 
almost flat within the patches. This flatness of the po- 
tential energy as bond-angle changes is the so called "de- 
generate valency" [111 [H[ , and two consequences follow: 
(i) the floppy modes of the central-force kagome lattice 
remain zero potential energy in the kagome lattice com- 
posed of tri-block patchy particles, and (ii) the kagome 
lattice is isoencrgetic with other lattices which have the 
same number of contacts per particle in their attractive 
patches. This is shown in Fig. [5J In addition, we shall 
call the area of the particle surface that is not coated the 
"non-attractive" patch. 

Surprisingly, this flatness in potential energy is lifted 
by thermal fluctuations via the so called "order by dis- 
order" effect [37l - l4"l| ] . A stable kagome lattice has been 
observed from the self-assembly of tri-block Janus par- 
ticles at finite temperature in the experiment discussed 
in Ref. [2(|. In this Paper we discuss a theory describ- 
ing this effect of how open lattices can be stabilized by 
entropy. 

In Sec. [IT] we set up the equilibrium statistical mechan- 
ics for lattices assembled from patchy particles and dis- 
cuss small fluctuations around the stable reference state. 
In Sec IIIII we discuss the effects of rotational and vibra- 
tional entropy of the patchy particles and how they con- 
tribute to the stability of open lattices. In Sec [IV] we 
summarize the results of our theory applied on the sys- 
tem of triblock Janus particles. In Sec [V] we discuss the 
agreement of our theory with experiment, possible im- 
provements, relations and applications to other systems. 
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FIG. 2. Examples of possible structures formed by triblock 
Janus particle of the same potential energy: (a) kagome lat- 
tice having bond angle 120°, (b) twisted kagome lattice with 
bond angle deviating from 120° (the twisted kagome and the 
kagome lattice is related by a floppy mode), (c) roman mosaic 
lattice consists of two types of particles with bond angles 120° 
and 150° (90° on the opposing side), and (d) hexagonal lat- 
tice with bond angle 120° but also contacts in non-attractive 
patch (red dashed lines) . The triangle denoted by orange dou- 
ble lines in (c) marks a corner-sharing-triangle which defines 
a basic unit in the structure. 



II. FREE ENERGY OF OPEN LATTICES 
ASSEMBLED FROM PATCHY PARTICLES: 
GENERAL FORMULATION 

To model the equilibrium statistical mechanics of a 
colloidal open lattice, we start from the generic form of 
partition function applicable to anisotropic particles, 

Z = e "bt [[dr-jdhj, (2.1) 
J j 

where H is the Hamiltonian, , hi are the position and 
orientation of the particle i, ks is the Boltzmann's con- 
stant, and T is the temperature. For a three-dimensional 
particle, the orientation n, represent the information of 
the polar and azimuth angle of the particle. If the parti- 
cle does not have symmetry for rotations around its own 
axes hi an additional variable is needed to fully describe 
its orientation. For simplicity we just consider particles 
with this symmetry and thus hi is sufficient. 

Now we specialize our discussion to the case spheri- 
cal shaped colloidal particles with inhomogeneous surface 




FIG. 3. Pair-wise interaction between two patchy particles. 
The green line denotes the vector connecting the centers of the 
two particles, tj — r,. The unit vector e»j = (ij — r<)/|rj — r,:|. 
The red arrows denote the orientations of the two particles, 
hi and hj. The surface elements crossed by the green line 
on the two particles dominate the interaction in Eq. 1)2.21) at 
short range. 

coatings (patches). The triblock Janus particles used in 
the experiment in Ref. (26j is one example of such par- 
ticles. For these particles the Hamiltonian of the system 
consists pair- wise potential energy of the form 

V(ri,hi,Tj,hj) = J dClidQjv(hi,hj,mi,mj, 

\(r t + am l /2) - (r,- + am,-/2)|), 

(2.2) 

which comes from the interactions v between the surface 
elements of the two particles, integrated over the solid 
angles of both of them. In this expression a is the diame- 
ter of the particle, rhi and rhj are unit vectors pointing in 
the directions of the solid angles Hi and Slj respectively, 
and thus | (r, + ami/2) — (rj + arhj/2) | is the distance be- 
tween the two surface elements. Figure [3] illustrate such 
a pair-wise interaction. If the range of interactions be- 
tween the surface elements of the two particles is much 
shorter than the particle size, and the two particles are 
sufficiently close, the surface elements crossed by the vec- 
tor pointing from the center of particle i to j dominates 
the integral and we can write 

V(l*i ; hi , Yj , hj ) —V i^Tli , 71 j , Cij , 6ij ; , 

|(r; + aey/2) - {v 3 - ae 4J /2)|) , (2.3) 
where is the unit vector pointing from particle i to j. 

We further specialize to patchy particles that exhibit 
short-ranged attractions when they face each other in 
their attractive patches and hard-sphere repulsion other- 
wise, and consider the free energy of equilibrium colloidal 
lattices assembled from these particles. For these patchy 
particles it is often practical to assume that (i) the depth 
of the attractive well is sufficiently greater than ther- 
mal energy (in the case of Ref. [26| the depth is approx- 
imately 7kgT), and (ii) the boundaries of the attractive 
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patch is sharply defined, i.e., the interactions are attrac- 
tive when two particles are facing each other in their 
attractive patch, and rapidly changes into hard-sphere 
repulsion otherwise. 

Thus for our consideration of the free energy of equilib- 
rium colloidal lattices, for which it is sufficient to consider 
only small fluctuations away from the stable state, we can 
make further simplifications of the problem. For a given 
lattice structure, we can identify all contacts between 
neighboring particles in their attractive patch, which we 
shall call "attractive bonds" , and in the calculation of the 
partition function only include those states in which these 
attractive bonds are not broken. All states with broken 
attractive bonds have a very small Boltzmann factor and 
can thus be ignored. For these attractive bonds, the po- 
tential energy can then be written into the form 

V(ri,ni,rj,hj) = v a (\ri ~ Tj\)w(eij,fk)w{—eij,hj) 

(2.4) 

where the factors 

!f if &ij passes through the 
attractive patch of particle i (2.5) 
oo otherwise 

and the central- force part of the interaction, v a (\ri —Tj\) 
is characterized by a hard-core repulsion, a short range 
attraction, and vanishing interaction when the distance 
between the centers of the particles is greater than the 
radius plus the range of attraction. 

There can also be contacts in the non-attractive sur- 
faces of the particles in the stable state (we shall call 
them "non-attractive bonds" ) , as we shall discuss in the 
case of the hexagonal lattice. These contacts are simply 
described by isotropic hard-sphere repulsions, given the 
fact that the stable state we start from is a local mini- 
mum in potential energy and these contacts will not turn 
into an attractive contact unless other attractive bonds 
are broken. Thus the potential of these non-attractive 
bonds is simply 

V(Ti,ni,Tj,hj) = v r (\ri - Tj\), (2.6) 

where v r describes a very short-ranged repulsion. 

With the form of potential energy (j2.4[) for the attrac- 
tive bonds, we can decouple the orientational degrees of 
freedom of different particles as follows. Because w — > oo 
if any of the particle has an orientation that one or more 
attractive bond is broken, and thus the corresponding 
Boltzmann factor vanishes in the summation of the par- 
tition function. We have 

J j 
= Je *5* ][*,({r,}./i,)./iy/H, (2-7) 



where the sum is over nearest bonds (mn) . The central- 
force Hamiltonian is given by 

i?CF({r l }) = ^ v a (\r m -r n \)+ ^ v r (\r m -r n \), 

(mn) a (mn) r 

(2.8) 

where (mn) a represent attractive bonds and (mn) r rep- 
resent non-attractive bonds. The factor ^ ({rj}, hj) = I 
if all of the attractive bonds of particle j remain in its 
attractive patch given its orientation, and if any attrac- 
tive bonds of j is broken. A two-dimensional example of 
Vfj function is shown in Fig.|4] It is clear that besides the 
positions of particles, this factor only depends on the 
orientation of one particle, regardless of the orientation 
of all other particles. Our assumption of small fluctua- 
tions near the equilibrium state allowed this decoupling. 
Therefore, from Eq. (|2 . T[) we can integrate out the orien- 
tational degrees of freedom of each particles individually. 
The result from each of these integration 

= / *3(to}.ftj)<*ni, (2.9) 

is proportional to the number of microscopic orienta- 
tional states allowed for particle j given the positions of 
this particle and its neighbors. It is appropriate to use the 
concept of microscopic orientational states here, because 
with fixed positions, the allowed rotation of one particle 
can be described by micro-canonical ensemble, since all 
allowed orientations are of the same energy. The result 
of this integral can be written in terms of the rotational 
entropy 

s ] =k B \nn j , (2.10) 

which adds up with the central-force part of the potential 
energy and lead to an effective Hamiltonian that depends 
on particle positions only 

-^Effective = #CF — T ^ Sj. (2-H) 
3 

The partition function is then 

C ^Effective 

Z= e fe s r [[drj. (2.12) 

J 3 

We shall discuss this calculation of the rotational entropy 
in detail for the case of triblock Janus particles below. 

III. EFFECTS OF ENTROPY: THE EXAMPLE 
OF TRIBLOCK JANUS PARTICLE 

In this section, we study the effects of entropy using 
the formulation discussed in Sec. HTJ We use the triblock 
Janus particles system as an example to show how en- 
tropy provides mechanical stability and induces differ- 
ences in the free energy of open and close-packed lattices 
which have the same potential energy. 
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(a) (b) 



FIG. 4. (a) An example configuration that gives 
$fj hj) = 1. In this configuration all bonds (green solid 

lines) connecting particle j go through its attractive patch 
(black area). The orange double arcs represent the allowed 
microscopic orientational states of particle j given the bonds, 
(b) An example configuration that gives ({fi }, n,) = 0. 
One bond goes through the non-attractive patch. The black 
dashed lines mark the boundary of the attractive patches. 
Rules for three dimensional configurations follow similarly. 

We will discuss two different type of triblock Janus par- 
ticles, corresponding to two recent experiments. The first 
type of triblock Janus particles have attractive patches 
that are elongated in the plane of the lattice [26[ . As a re- 
sult for the lattices we discuss in this Paper, the kagomc 
and the hexagonal lattices, in which each particle has 
4 attractive bonds, the orientational fluctuations of the 
particle is confined to two dimensions in the calculation 
of partition function (|2.7[) . The second type of triblock 
Janus particles have circular shaped attractive patches, 
and the allowed orientational fluctuations are three di- 
mensional [13]. We will discuss both cases in the follow- 
ing. 

A. Rotational entropy 

1. Triblock Janus particles with elongated patches 

These triblock particles can be considered as two- 
dimensional disks that assemble into two-dimensional lat- 
tices, and their orientations are characterized by only the 
polar angle 9i. 

As shown in Fig. [5l to keep all attractive bonds of par- 
ticle i in its attractive patch, its orientation 9i is confined 
by the bond angles, a, a', 71, 72 (here and in the follow- 
ing discussion in this subsection we drop for subscript i 
for convenience), between the 4 attractive bonds. It is 
clear that there are only 3 independent variables because 

a + a' + 71 + 72 = 2n. (3.1) 

The values of these angles in the stable state are 

ot = a=—, 71=72=3, (3.2) 
and they deviate from these values at finite temperature. 




FIG. 5. A triblock Janus particle with rotations confined in 
two-dimensional space, due to the elongated patch shape (not 
shown in this figure). The patch size cj> and the bond angles 
a, a', 71, 72 are also shown. The 4 green solid lines represent 
the 4 attractive patches of the particle, and the red arrow 
represent the orientation of the particle. The black dashed 
lines represent the boundaries of the attractive patches. 

As we found by comparing to experimental data, in lat- 
tices assembled from triblock Janus particles, the stiffness 
against bond-angle fluctuations is much smaller than the 
stiffness against bond-length fluctuations 34]. Because 
the two pairs of nearest neighbors of the 4 attractive 
bonds one triblock Janus particle (in Fig. (SJ the top pair 
and the bottom pair) are themselves nearest-neighbor 
pairs, it is straightforward to realize that the fluctuations 
in the angles 71 and 72 are much smaller than those of 
the angles a and a' . Therefore, it is reasonable to make 
the simplifying assumption of ignoring the fluctuations 
in 71 and 72, so that 

7T 

71=72= (3.3) 

Then we have 

4-7T 

a + a' = —, (3.4) 

and there is only one bond-angle variable a. We also 
show the calculation of rotational entropy with fluctu- 
ations in 71,72 allowed in Sec. IIII A 31 The calculation 
is considerably more complicated but the result is very 
close of that of the current simplified model. 

Given a we consider the allowed orientations 6 of the 
particle. Choosing the polar coordinate so that at 8 = 
the "north pole" of the particle is aligned with the axes of 
reflection symmetry of the 4 attractive bonds, as shown 
in Fig. [5j The total allowed rotational microscopic states 
of the particle, as introduced in Eq. (|2.9p . is thus propor- 
tional to 

_£„e(. + ,-^) 

x6(^-(0 + n-d>)) (3.5) 
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FIG. 6. 3D rotation of the sphere. Green lines denote the 
bonds, and the red arrow denote the north pole of the particle. 



where <p is the angle characterizing the half patch size 
as shown in Fig. [5J and without loss of generality we 
assumed a < a'. We used the Hcavisidc step function, 
6, to enforce our condition of the >F functions, that the 
allowed orientations 9 need to keep all attractive bonds 
in the attractive patches of this particle. Introducing the 
bond angle fluctuation variable 



we have 



2tt 
~3~' 



(3.6) 



(3.7) 



where the absolute value operation on Act comes from 
releasing the assumption that a < a' . This result leads 
to the rotational entropy of the particle given bond angle 
change A a , 



s = ks In 



|Aa| 



(3.8) 



This rotational entropy is shown together with that of 
the particles of circular patches in Fig. [7] For any given 
set of particle positions {r^} all the bond angles are given 
and we can calculate the rotational entropy Si for each 
of the particles. 



2. Triblock Janus particles with circular patches 

In contrast to our discussion on Janus particles with 
elongated patches, particles with circular patches need 
to be described using three-dimensional rotations with 
polar angle 9 and azimuth angle tp. 

For small fluctuations around the stable state, we still 
require that the 4 attractive bonds remain in the attrac- 
tive patch. Given that the 4 attractive bonds are in the 
same plane (for two-dimensional kagome and hexagonal 
lattices), and 71 = 72 = 7r/3, we can write down the unit 
vectors pointing in the directions of these 4 attractive 



bonds 
Mr = 
M 2 -- 
M 3 = 

Mi = 



{sin(f-f),0,cos(|-f)}, 
r /7T a\ /7r a\ ~i 



7r a 

— sm 1 

2 2 

7T a' 

— sm [ 

2 2 



,0,- 



7T 
' 2 



a' 






) 




a' 


+ 


2 



(3-9) 



where a and a' are the bond angles on the left and the 
right and they are related through 



4tt/3. 



(3.10) 



We then keep these 4 bonds fixed in space and let the 
particle rotate around its center, and calculate the rota- 
tional microscopic states of this particle given that all 4 
attractive bonds are kept in its attractive patches. The 
direction of the "north pole" of the particle 



N = {sin 9 cos ifj, sin 9 sin ip, cos 9} 



(3.11) 



fully describes the orientation of this particle, because 
the patches of the triblock Janus particles are symmetric 
with respect to rotations around its north pole. If the 
patches are of more complicated patterns, more variables 
will be needed to describe it. 

Again, without losing any generality we assume that 
a < a' . The condition that all four bonds remain in the 
attractive patch is then 



N ■ Mi > cosi 
-N ■ M 2 > cos ( 



which simplify into 



(X OL 

cos 9 sin — ± sin 9 cos ib cos — > cos < 
2 2 ~ 



(3.12) 



(3.13) 



The total number of rotational microscopic states is thus 
proportional to 



d6 
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a a 

x cos 9 sin h sin 9 cos ib cos — 

2 2 
a ol 
x ( cos 9 sin — — sin 9 cos ip cos — 



COS( 



(3.14) 



Equation (|3.13j) indicates two critical value for 



a , 7T 

2 +4> -r 



arccos 



cos ( 



(3.15) 



that for < 9 < 9\ the equality (|3.13[) is satisfied for all 
values of %p, and for 9\ < 9 < 9 2 the equality (|3.13[) is 
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only satisfied for some values of tp, in this case, around 
%jj = ±7r/2. Therefore the number of rotational states in 
the three-dimensional rotation case is proportional to 



(3.16) 



with 



= 2tt(1 - cos6»i) 
(2) = / 



L2 



— arccos 



sin 6 cos j 



(3.17) 



We calculated this number of rotational states and plot 
the results in Fig. [7] in the paper. It has no qualitative 
difference from the 2D case. 




FIG. 7. Rotational entropy of the triblock Janus particles 
with elongated patch (and thus 2D rotations), denoted by 
solid lines, and circular patch (and thus 3-D rotations), de- 
noted by dashed lines, as a function of the bond angle differ- 
ence Aa at different patch size <j> as listed in the legend. The 
curves from outer to inner represent the change of rotational 
entropy s — so in cases of <j> = 55° , 45° , 35° respectively, where 
so denote the rotational entropy when Aa = 0. 



Correction to rotational entropy allowing bond length 
fluctuations 



A more realistic description of the rotational entropy 
should allow the fluctuations in bond length. In fact in 
Sec. IIIIBI we shall discuss the contribution of these vi- 
brational modes to the free energy. These bond length 
fluctuations lead to fluctuations in the angles 71 and 72, 
which we took to be fixed value 7r/3 in our previous dis- 
cussions. For convenience we define the deviations of the 
angles 



a = 27r/3 + Aa, a' = 2tt/3 + Aa', 
71 =7r/3 + A71, 72 = 7t/3 + A7 2 , 



(3.18) 



and they satisfy the equation 



Aa + Aa' + A71 + A72 = 0. 



(3.19) 



Thus we just choose Aa, A71 , A72 to be the independent 
variables. We now calculate the rotational entropy in 
the two-dimensional case. The case of three-dimensional 
rotations follow similarly but involve more complicated 
derivations. 

Similarly, to calculate the number of rotational micro- 
scopic states, we that all the 4 bonds stay in the at- 
tractive patch. We arrive at the following form for the 
rotational entropy 




2{<f> - 7r/6) + Aa 
2(0-tt/6)- A72 
2(cf> - tt/6) - A 7l 
2(0 - tt/6) - A72 - A 7l 



if Aa + A71 < and Aa + A72 < 0, 
if Aa + A71 < and Aa + A72 > 0, 
if Aa + A71 > and Aa + A72 < 0, 
Aa if Aa + A71 > and Aa + A72 > 0, 



for the range of — 2(0 — tt/6) < Aa < 2(0— tt/6) — 71 — 72 
(the number of rotational states vanishes for Aa outside 
this range). It can be easily understood as follows. In the 
case of Aa + A7! < and Aa + A72 < 0, it is the two 
bonds on the right that determine the allowed rotations. 
In the case of Aa + A71 < and Aa + A72 > 0, it is 
the two bonds on the bottom that determine the allowed 
rotations. In the case of Aa+A7i > and Aa+A72 < 0, 
it is the two bonds on the top that determine the allowed 
rotations. In the case of Aa+A7i > and Aa+A7 2 > 0, 
it is the two bonds on the left that determine the allowed 
rotations. These different cases are depicted in Fig. [5] 



The deviations 71 and 72 are limited by the central- 
force interactions between nearest neighbors and are thus 
small. We can then assume that I71I, I72I -C 2(0 — 7r/6) 
for the range of parameters we are interested in. 



B. Vibrational entropy 

In this section we discuss the next step: integrating out 
the positional degrees of freedom in the partition func- 
tion and arrive at the free energy for given lattices. We 
shall find that although open and close-packed lattices 
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(b) 




(c) 



(d) 



FIG. 8. Bond angle configurations at given 71 and 72 in the 4 
situations in Eq. (pT20)l . (a) AQ + A71 < and Aa + A7 2 < 0. 
(b) Aa + A71 < and Aa + A72 > 0. (c) Aa + A71 > and 
Aa + A72 < 0. (d) Aa + A71 > and Aa + A72 > 0. 



can have the same potential energy, due to their differ- 
ence in vibrational entropy, their free energy is different. 
Open lattices, which exhibit lower free energy, can be 
selected in the system of patchy particles due to these 
entropic effects. 



1. Effective Hamiltonian, mechanical stability, and 
harmonic approximation 

The rotational entropy discussed above can be plugged 
into the equation for the effective Hamiltonian (|2.11l) . 
which is a version of the Ginzburg-Landau free energy 
with orientational degrees of freedom integrated out but 
positional degrees of freedom of the particles kept. The 
rotational entropy is related to the orientational micro- 
scopic states we calculated above through 



= kn hi O 



3 ' 



(3.21) 



where VI j is the microscopic states calculated for the j-th 
particle following the methods we discussed. 

The resulting -/^Effective is an effective Hamiltonian that 
only depends on particle positions. In particular, with 
these rotational entropy terms, more constraints (in the 
Language of Maxwell's counting) are being introduced 
into the system comparing to the original Hamiltonian 
which has potential energy only. For example, in the case 
of kagome lattices assembled from triblock Janus parti- 
cles, the rotational entropy terms introduce constraints 
on the bond angles, which were the floppy modes in the 
potential-energy only Hamiltonian, as we discussed be- 
fore. Thus as a result the rotational entropy may in 
general provide mechanical stability to the open lattices, 



which are not stable at T = from potential energy con- 
siderations. 

To obtain the free energy for each lattice structure, 
we need to also integrate out the positional fluctuations. 
In the mean time, it is also important to characterize 
the vibrational modes in these assembled lattices, which 
provide a valuable information in studying the mechan- 
ical response of the assembled structure. In order to 
achieve this, we choose to use the harmonic approxima- 
tion, i.e., keeping the effective Hamiltonian to quadratic 
order in displacement vectors Uj of particles, which facil- 
itate analytic calculations and is convenient in extract- 
ing asymptotic behaviors. A full calculation of the free 
energy involves more realistic modeling of the interac- 
tion potentials, including the hydrophobic and screened 
electrostatic interactions between the particles, which are 
anharmonic. We shall investigate this full calculation nu- 
merically in our future study. 

In order to obtain a harmonic form of the effective 
Hamiltonian, we shall make approximations on the ro- 
tational entropy so that it can be written into a form 
involving only quadratic terms in vij. Given the special 
nature of the patchy particles, as shown in Fig. the 
rotational entropy is linear at small Aa but nonlinearity 
rises at slightly larger Aa. In practice the fluctuation of 
the bond-angle a can be relatively large owning to the 
low stiffness of the corresponding vibrational mode, as 
we shall discuss later. Thus, it is a good approximation 
to require that the second moment ((Aa) 2 ) of the bond- 
angle fluctuations controlled by the rotational entropy in 
the original form should be kept the same in the harmonic 
approximation. For the case of the elongated patches we 
have 



((Aa) 



J(Aa) 2 e s / kB dAa 
J e s / kB dAa 

§(*-r 



(3.22) 



where s is given by Eq. (|3.8j) . This average can either 
be viewed as weighted by the effective Gibbs factor from 
the entropic term of the rotation of the given particle 
in the effective Hamiltonian (|2.11|) . or as weighted by 
the number of microscopic states of the particle, because 
e s/k B — fi The corresponding harmonic form with the 
same second moment is then 



-Ts 



3k B T 



(3.23) 



where the resulting bond-angle stiffness k is proportional 
to ksT because it is of entropic origin. For the case of 
circular patches we can do a similar calculation, the result 
will be a modified value of the bond-angle rigidity k. 
The case discussed in Sec. IIII A~3l is different, because of 
the additional dependence on the bond angles 71 and 72. 
We discuss the harmonic approximation of the rotational 
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entropy of this case in App. \K\ and the result is 

-Ts^-^Aa + ^±^l)\ (3.24) 
where n is the same as given in Eq. (|3.23[) . 



2. Dynamical matrix of the lattices 

With the harmonic forms of the rotational entropy de- 
rived above, in order to write the whole the effective 
Hamiltonian (|2.11[) into a harmonic form and thus ex- 
press it in terms of a dynamical matrix, we also need to 
characterize the central force interactions i?cF and make 
harmonic approximations on this part too. 

The central force interactions between particles facing 
each other in their attractive patches are short-range at- 
tractions. In the experiment of Refs. [26|, |34| they are of 
hydrophobic nature. In order to make a harmonic ap- 
proximation we expand this short-range attraction near 
its minimum to quadratic order, which is equivalent to a 
harmonic spring, and the potential can be written as 



Rj — Ri\ — I 



(3.25) 



where we assume the particles undergo a displacement 
from original positions to displaced positions Ri 



-> Ri = n 



(3.26) 



and k is the spring constant, I is the rest length. 

For particles not facing each other both in their at- 
tractive patch the central force interaction is essentially 
a hard sphere repulsion. In practice the repulsion is of a 
very small range. In the experiment of Refs. (2(| [H[ this 
repulsion is of the nature of strongly-screened electro- 
static repulsion with the range of repulsion comparable 
to surface roughness. We can also expand this interaction 
around the stable inter-particle distances in the lattice to 
quadratic order and approximate it as a harmonic spring. 
Because this interaction is purely repulsive in nature, the 
expansion in terms of distance variations exhibit a linear 
term, which contribute to the internal stress of the lat- 
tice but do not change our analysis of vibrational modes. 
Similar to the attractive bonds, the potential of these 
non- attractive bonds can be written as 



V, 



2 



Rj — Ri | — I 



(3.27) 



where k r is the spring constant for the non-attractive 
bonds, I is the rest length. 

We can then expand the central-force potentials in 
Eqs. (|3.25l3.27[) in the displacement vectors Ui and keep 
to quadratic order. The same procedure needs to be ap- 
plied to the rotational entropy terms through the relation 
of bond angle with displacement vectors in 2D 



Aa, ; 



-(fhi x u h i - hj x Uij) ■ z 
a 



(3.28) 



where h,i,j labels the particles forming the bond angle, 
fhi and f^i are the unit vectors pointing from h to i, and 
from i to j in the undeformed state. The unit vector 
z points in the third dimension out of the plane. The 
diameter of the particle, which is the same as the distance 
between nearest neighbors, is denoted by a. The same 
formula applies to A71, A72 too. 

Thus the effective Hamiltonian, which contains both 
the central-force terms and the rotational-entropy terms 



k 



Effective 



J2 -(i^-^i-z 



IE (A 



^oti + g- 



A7 M + A72,. 



(3.29) 



where the summation (ij),a is over attractive bonds, 
(ij),r is over non-attractive bonds, i in the third term 
is over all particles. The parameter g = if we adopt the 
simple form (|3.23|) of the rotational entropy and g = 1 if 
we adopt the form (|3.24|) which allows bond-length fluc- 
tuations. In calculating the figures in this Section and 
the next, we used 5 = 1. The result will be qualitatively 
the same if we use g = 0. This effective Hamiltonian 
of the lattice can be readily written in the formalism of 
dynamical matrix as 



-^Effective = ~ ^2 UtDt^tUy, 



(3.30) 



where I labels unit cells each containing m particles, and 
Ui denote the 2m-dimcnsional displacement vector of the 
particles in the unit cell in two dimensions. In App. IB II 
we derive the dynamical matrix for the hexagonal and 
the kagome lattices. 



3. Normal modes of the lattices 

Normal modes of the two lattices can be directly de- 
rived from their dynamical matrices D^^ as cigenmodes 
of the matrices. For periodic lattices this analysis is most 
conveniently done in momentum space using D,^ . 

For the hexagonal lattice, because D,^ is 2 x 2, there 
are two branches of modes as shown in Fig. [9^, corre- 
sponding to the two acoustic phonons in two dimensions. 
The frequencies of both of the two branches are controlled 
by fc, k r , and k. In this plot we used the values of fc, 
k r , and k from fitting to the experiment as discussed in 
Ref. [34| and also the case with k = 0. It is straightfor- 
ward to see that the difference in the frequency in the 
case with k = and k > (but small compare to ka 2 ) is 
small and only quantitative in the hexagonal lattice. 

For the kagome lattice, there are six branches of modes 
as shown in Fig. ^jp, because each unit cell consists 
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tional to s/k. Thus the two cases of k = and k > are 
qualitatively different for the kagome lattice. In partic- 
ular, because of the lnw contribution to the free energy, 
as k becomes small, the free energy of the kagome lattics 
is significantly lowered, as we will discuss in Sec. IIVI 

In addition, the eigenvalues (square of the frequencies) 
of the six modes of the kagome lattice at q = are given 
by 



12k 
{0,0,^-,3fc 



3(1 + 3g 2 )n 
2d? ; 



3fc 



3(1 + 3g 2 )n 
2^2 : 



Qk} 
(3.31) 



(a) 



where the floppy mode is the third one in the list. 





FIG. 9. Phonon dispersion relations for the hexagonal lattice 
(a) and the kagome lattice (b). The dispersion relations are 
plotted in the first Brillouin zone. The black solid curves rep- 
resent the case with k = 0, k = 2300ksTa~ 2 and k r = 20fc. 
The red dotted curves represent the case with « = 33fcsT, 
k = 2300fcsTa -2 and k r — 20k. The hexagonal lattice show 
relatively high frequency in some part of the dispersion rela- 
tion plot because of the high value of k r . The first Brillouin 
zone for the hexagonal and the kagome lattices are shown in 
(c). 



three particles. As discussed in Rcf. |43j, there are two 
acoustic branches and four optical branches of phonons. 
In the case of K = the lowest optical branch drops 
to zero frequency (upon hybridization with transverse 
phonons) along the TM direction in the first Brillouin 
zone. This branch corresponds to the floppy modes of the 
lattice, which are of zero-frequency if there is only central 
force interactions between nearest neighbors, as shown in 
Fig. [TJ For k > this branch gains a gap that is propor- 



IV. RESULTS: COMPETITION BETWEEN 
OPEN AND CLOSE PACKED COLLOIDAL 
LATTICES 

A. Free energy difference 

At finite temperature and in thermal equilibrium, the 
relative stability of different lattice phases is determined 
by their free energies per particle /. The free energy 
per particle of a lattice can be directly evaluated using 
the partition function (|2.12p and the effective Hamilto- 
nian p.30[) we developed above, 



/ = -{Nmy 1 k B T\nZ 

k B T , „ „ d, m , ( 2ixk B T 

= 77^— m DetD - -k B T\n ^~ 

2Nm 2 V a 2 



(4.1) 



where the second term is a constant which we shall denote 
as C. The determinant is taken over the mNd dimensions 
of the tensor -D"^, where label the unit cells and a/3 
label the Cartesian indices. 

We can further write this free energy using the 
dynamical matrix in momentum space as defined in 
Eqs. (|B4IB5IB6I) . in which it is diagonal in q and con- 
venient to evaluate 



/ = 



In 



k B T 
2Nm 
k B T 
2Nm 
k B T 

2Nm 
k B T v 

2m V 



In DetD + C 



(Vv Q ) 



-Nm 



DetD 



C 



In Det£> + C 



^lndctD 9 -+C', 



(4.2) 



where the determinant in the last line is only taken in 
the 2m dimensional space at a given q. In the continuum 
limit we can write this into 



/ 



k B T 
2m 



IBZ 



f-^-lndetD„- 

(27T V 



c 



(4.3) 



where the integral is over the first Brillouin zone. In what 
follows we specialize to the free energies of kagome and 
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FIG. 10. The increase of kagome lattice free energy per par- 
ticle as a function of the bending stiffness k. Blue dots repre- 
sent numerical integrals and the red line represent fitted form 
(/k - /xU=o) /(k B T) ~ b y/K/(ka 2 ) with bo = 2.55. 



hexagonal lattices assembled from triblock Janus parti- 
cles. 

Firstly, we consider the kagome lattice. In the n = 
case, the floppy modes along TM lines lead to a logarith- 
mic divergence in the lndetS^ terms in the integrand. 
However the free energy which is after the integral over 
the first Brillouin Zone does converge to a finite value, 



/if U=o = k B T 



- In k + - In - 
2 2 V 4 



C, 



(4.4) 



where Jk denotes the free energy per particle of the 
kagome lattice. It is noteworthy that the dependence 
on k is through an additive term, because in the case of 
k = all terms in 2) is proportional to k and thus k can 
thus be factorized out. This term can be combined with 
C leading to a term baring the dimension of energy 



2 B \2nk B T 



so that 



Ik U=o 



2 \4 



C. 



(4.5) 



(4.6) 



For the comparison between different lattices consisting 
the same type of patchy particles, this term is a lattice- 
independent constant and will cancel out in the free en- 
ergy differences. 

For the case of k > but small, the lead ing order 
correction to fx is proportional to y/n/ (ka 2 ), as shown 
in Fig. 1101 This non-analytical dependence on k comes 
from the logarithmic divergence in In a;, which originates 
from the mechanical instability of the kagome lattice. 

Secondly, we consider the free energy of the hexagonal 
lattice. In the case of k = 0, the lattice is stable and 
there is no divergence in the integrand at all. Owing to 
the extra parameter k r associated with the non-attractive 
bonds, the free energy per particle of the hexagonal lat- 
tice //fU=o depends on both k and k r . In particular, 
given the logarithm, k can be factorized out contributing 
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FIG. 11. (a) The difference between the hexagonal lattice and 
the kagome lattice free energy per particle as a function of 
k r /k at k = 0. At k r /k — and k — there is no difference 
between these two lattices so the curve goes to zero. For 
kr/k of order 1 or greater, the difference between the two free 
energies is of the order of IksT. (b) The hexagonal lattice 
free energy per particle as a function of «/ (ka 2 ) at k r /k — 20. 
It is clear that in this case the correction due to k is linear. 



to the same constant C as in the kagome lattice. The 
other term, instead of being a constant, bares a depen- 
dence on the ratio k r /k. In the limit of k r — this just 
becomes the same problem as the kagome lattice with 
K = 0. In real systems typically k r is of the same order 
of magnitude as or g reater strength than k. In the ex- 
periment in Ref. [34[ it is found that k r /k ~ 20. The 
dependence of /frU=o on k r /k is shown in Fig. [TTTa). 

For the case of k > 0, fjj increases with k with a 
leading order correction proportional to n/(ka 2 ) instead 
of y/ kJ 1 (ka 2 ) (except for the case in which k r is very 
small), because the hexagonal lattice is not affected by 
the marginal stability as is the kagome lattice. This is 
shown in Fig. [TTT bl . 

Therefore, in calculating the difference between the 
free energy per particle of the hexagonal and the kagome 
lattices, the constant term C cancels out, and the re- 
sulting difference depends on the two ratios n/(ka 2 ) and 
k r /k. In particular, A/ = fjj — fx decreases with 
n/(ka 2 ), and increases with k r /k, as shown in Fig. 1121 
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FIG. 12. Free energy difference between the hexagonal and 
the kagome lattices as a function of K,/(ka 2 ) at k r = 1 (lower 
blue curve) and k r = 20 (upper red curve). 
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FIG. 13. Equilibrium phase diagram showing the boundary 
between the hexagonal and the kagome lattice determined by 
the equal Gibbs free energy per particle line. 



From this calculation we see that the free energy dif- 
ference is significant, i.e., of order fcgT, in two situa- 
tions: (i) n/{ka 2 ) <C 1, or (ii) k r /k ^> 1. This provide us 
with the criterion of the selection of the kagome lattice 
over the hexagonal lattice by free energy considerations. 
For the case of k r /k not significantly greater than unity, 
the selection of the kagome lattice is only effective when 
n/(ka 2 ) ^ 1. Because the bending stiffness k originates 
from the rotational entropy as shown in Eq. (|3.23j) . this 
condition leads to a characteristic patch size 



7T 

G 



3k B T 
2ka 2 



(4.7) 



above which the selection of the kagome lattice is effec- 
tive (if k r /k not already significantly greater than unity). 
This reveals the perhaps surprising effect that the selec- 
tion of open lattices can be enhanced by having exces- 
sively large attractive patch (34|. 



B. Phase diagram 

At positive pressure there exist a first-order transition 
from the kagome to the hexagonal lattice, because the 
hexagonal lattice is of higher packing fraction [34| . This 
transition can be described by transform from the statis- 
tical ensemble of fixed lattice structure into the ensemble 
of fixed pressure through a Legendre transformation 



9 = f + pv 



(4.8) 



where g is the Gibbs free energy per particle, p is the 
pressure, and v = V/N is the mean area occupied by 
one particle in the given lattice. At fixed p in equilib- 
rium, the phase with lower g is more stable. Therefore 
we can determine the boundary between the kagome and 
the hexagonal phases using the equal g line. 

For a qualitative description we make further simplifi- 
cation that within a phase the change in the lattice con- 
stant as a function of pressure can be ignored, thus the 



only volume change comes from the difference between 
the two lattice structures. In particular the value of v in 
the kagome and the hexagonal lattices are respectively 



v K 



2a? 



V H 



(4.9) 



and we define Av = Vjj — vx- 

We also make the simplifying assumption that the stiff- 
ness parameters k, k r and k also remain the same as pres- 
sure changes, because the lattice constant is unchanged. 
Thus the difference in Hclmholtz free energy between the 
two lattices, A/ = fu — fx, is the same at different pres- 
sures. In addition, the potential energy is the same in the 
two phases, because each particle have 4 attractive bonds 
in both phases, and as a result, wc have A/ = — TA.s 
where As is the entropy difference per particle between 
the two lattices. 

Therefore the equal g line, which describes the phase 
boundary, correspond to the equation 



P 



T- 



Av' 



(4.10) 



This leads to the phase diagram shown in Fig. [T3] If this 
phase diagram is plotted in the traditional p — T plane, 
the phase boundary will be a straight line, of which the 
slope positively depend on the patch-size. In addition, at 
high temperature we expect that both lattices will melt 
into a fluid phase, which is not considered in this theory. 



V. DISCUSSION 

In this Paper we construct a generic theory based on 
equilibrium statistical mechanics and lattice dynamics for 
the self-assembly of periodic lattices from patchy par- 
ticles. We discuss the entropic effects in open lattices 
formed by patchy particles. In particular, wc show that 
the rotational entropy can provide mechanical stability to 
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open lattices and vibrational entropy can lower the free 
energy of the open lattice relative to that of close-packed 
lattices. These effects are essential to the stabilization 
and selection of open lattices which are not stable in pure 
potential energy considerations. 

The predictions of our theory agrees well with exper- 
iments on the triblock Janus particle system [34j]. In 
particular, the transition between the hexagonal and the 
kagome lattices predicted by this theory has been con- 
firmed in experiment with increased lateral pressure. The 
generic dependence of the phase boundary between the 
open and close-packed lattices on temperature and patch- 
size has yet to be confirmed experimentally. 

This stabilization and selection effect of entropy in 
patchy particle systems is an example of the "order by 
disorder" effect, in which degenerate states at T = gain 
different free energy due to thermal or quantum fluc- 
tuations [37M4T1 | . Related examples of this effect has 
been discussed in homogeneous colloids. For example, 
close-packed lattice structures can be favored in hard- 
sphere homogeneous-colloid systems by entropy, because 
the regular lattice structure allows more "rattle room" for 
the particles [H[. For patchy colloids, as discussed here, 
the cntropic effect is even more interesting owing to the 
nontrivial rotational degrees of freedom. 

In order to arrive at the form of the effective Hamil- 
tonian that is suitable for harmonic lattice dynamics, a 
series of approximations have been made in this theory, 
and the major ones are: (i) the pair-wise interactions is 
very short-ranged and is of step- function nature, i.e., flat 
in attractive and non-attractive surfaces areas and only 
changes abruptly as particle orientation is altered across 
the boundary; (ii) both the central-force inter-particle 
potential and the effective bending potential term from 
rotational entropy can be approximated as harmonic in 
displacement vectors u of the particles. 

The approximation (i) is made so that the orientational 
degrees of freedom of different particles can be decoupled. 
For the experiment in Refs. (2(| [34| the salt concentra- 
tion is high enough to screen the length-scale of electro- 
static repulsion down to the scale of surface roughness 
and the approximation (i) is valid in these cases. For ex- 
periments in which the pair-wise potential is not so flat 
as a function of particle orientations, the orientational 
degrees of freedom of different particles can not be de- 
coupled. Furthermore, the kagome lattice may directly 
gain some stability from potential energy. Nevertheless, 
typically for patchy particles this potential energy varia- 
tion as a function of bond angles is not strong enough to 
actually stabilize the lattice, and the effect of rotational 
entropy on the stability may still dominate. 

The approximation (ii) validates the application of lat- 
tice dynamics to the effective Hamiltonian of the system 
after orientational degrees of freedom are integrated out. 
We made this approximation to obtain qualitative predic- 
tions of the free energies of different lattices, which reveal 
the generic feature of the phase diagram without requir- 
ing complicated calculations. More accurate calculation 



of the free energies of the lattices can be done numer- 
ically with more realistic modeling of the inter-particle 
potentials. We do not expect qualitative difference in 
the conclusion by replacing the harmonic approximation 
with more realistic potentials. 

This theory is readily gcncralizable to more cases 
of self-assembly of patchy particles, including three- 
dimensional lattices. From the entropic stabilization 
mechanism discussed in this theory, it is expectable 
that in three dimension, transitions between close-packed 
face-centered-cubic and open pyrochlore/perovskite lat- 
tices occur in the tri-block patchy particle system, in a 
similar fashion as the hexagonal/kagome transition in 
two dimensions. This cntropic effect provides a mech- 
anism in which regular periodic lattices are automati- 
cally favored with remarkably simple designs of building 
blocks. 
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Appendix A: Harmonic approximation of the 
rotational entropy 

In this section we consider the Harmonic approxima- 
tion of the rotational entropy 

s = fc s lnfi, (Al) 

with fi given in Eq. (|3.20p . Due to the dependence on the 
angles 71 , 72 it is not possible to always keep the second 
moment ((Aa) 2 ) invariant in this approximation. There- 
fore we need to generalize our method of approximation. 

Here we show that the harmonic approximation we de- 
fined for the simple case of defined in Eq. (|3.7j) can 
also be obtained through the following more generic pro- 
cedure. We firstly define the Fourier transform of the 
statistical weight e s > kB 

y(q) = J e s ' kB e - iqAa dAa. (A2) 

Then we can take the logarithm of y{q) and keep to 
quadratic order in q. This leads to a Gaussian approxi- 
mation of y(q) as 

y(q) =y(0)e-£. (A3) 

The inverse Fourier transform of y(q) leads us to a Gaus- 
sian approximation of the original statistical weight 

e S/k B = f ^ly(^ q \ e i<lAa = ( A4 ) 
J 27T 
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It is straightforward to show that following such a pro- 
cedure of approximating a statistical distribution into 
a Gaussian distribution, the second moment ((Aa) 2 ) is 
kept invariant. The second moment can be calculated 
using the full Fourier transform y{q) as follows 



and 



((Aa) 2 ) = J dAa(Aa) 2 e s/kB / J dAc 



s/k E 



y"(Q) 
2/(0) ' 



(A5) 



Thus it is clear for this second moment we only need to 
keep to 0(q 2 ) on the exponent of y(q), which is y(q), and 
higher order terms will not contribute. 

Now we can use this method to calculate the harmonic 
approximation of Q given in Eq. (|3.20j) , which applies to 
the case with bond-length fluctuations allowed. We have 
the Gaussian form of the Fourier transform 



y(q) = y(o)e' 



/3^V+6i(A T1 +A 72 )q 



(A6) 



in which we also kept to quadratic order in A71 , A72 
because they are small fluctuations. Transforming back 
we arrive at 

rp K ( K , A71 + A 72 \ 2 

-Ts^--[Aa+ j , (A7) 

where k is the same as given in Eq. (|3.23|) . 



e\ = cos[(n — 2)7r/3]e x + sin[(n — 2)7r/3]e y 
= — sin[(n — 2)7r/3]e.,; + cos[(n — 2)7r/3]e. y 

(B3) 

are the unit vectors along or perpendicular to the bonds 
as shown in Fig.ll4h. The lattice constant of this hexago- 
nal lattice is equal to the bond length a. The momentum 
space dynamical matrix is related to the real space one 
through 

t,t' 

= ^2 S e-** 1 ***' D ?>a r (B4) 
<?,<?' 

following the convention of Fourier transform on lattices 
in Ref. [46|. The factor «o,h = V3a 2 /2 is the area of the 
unit cell for the hexagonal lattice. Using translational 
invariance the dynamical matrix in momentum space can 
be written into 



D <?\<r 



Vv Q . h Sf-^Df 



with 



(B5) 



(B6) 



Appendix B: Dynamical matrices and normal modes 
of the effective Hamiltonian 



1. Constructing the dynamical matrices 

a. The hexagonal lattice 



The bending part of this dynamical matrix is different 
from the version used in Ref. [45|, [47} for the reason that 
here each particle just has one bending energy term, and 
the two bonds forming this angle is 120° in the reference 
state rather than 180° as for the filaments. 



Starting from the effective Hamiltonian as given in 
Eq. (|3.29[) . we can construct the dynamical matrix of 
the hexagonal lattice following the method discussed in 
Ref. H|, 

^hexagonal =fc _ + g. ^s _^ + kJ^B^ 

+ ^r ?+ - 9 — 2 — 



where 




5' if =(l-e-^)4 

B\ = (1 - e"^ )ei + (e - **^ - l)e^ 
f* f = -(1 - e-^)e± + (e-^ s 3 - 1)^ 

= (1 - )ei - {e-^l - l)e£ (B2) 



(Bl) 



b. The kagome lattice 

The dynamical matrix for the kagome lattice can be 
constructed in a similar way. Because there are 3 parti- 
cles in one unit cell of the kagome lattice, the displace- 
ment vector Ug is 6 dimensional representing the 2d mo- 
tion of the 3 particles, 

U[ = {Ut,l,x, u Z,l,y, U£,2,x, u e,2,y, Ut,3 lX , Ui,3,y}- (B7) 

To denote the directions of the bonds in the lattice, we 
define the unit vectors 

e| = cos[(n + l)2n/3]e x + sin[(n + l)27r/3]e y 
e£ = - sin[(n + l)27r/3]eai + cos[(n + l)27r/3]e y 

(B8) 

with n = 1,2,3 for the unit vectors along and perpendic- 
ular to the bonds, as shown in Fig. [T4"b . 
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FIG. 14. The architecture and vectors in hexagonal (a) and 
kagome (b) lattices. The black numbers 1, 2, 3 marks the three 
unit vectors e| in the two lattices. The red disks with white 
numbers on them in (b) label the 3 particles in a unit cell in 
the kagome lattice. 



with 



and 



B 2,q 



1 1,9 ■ 

L 3,q 

f b - 

4,9 ' 

f b - 
1 5,q 

f b - 



Kq = {-ef,ef,0,0}, 
B s 24 = {0,0,-44}, 

B^= {4,0,0,-4}, 
Bl, f ={ef,-e 2 ^"ef,0,0}, 



^iq-e\A\ 



}■■ 



Bl ? ={-e 2 ^el0,0,e^, 



= {- 



5L- = {-^ 



e 3 > e 2 i 



'1 ' °3 



{e 1 , e 2 , e 2 }, 

{^3 i ^2 i — &2 ~ 63 }, 



= {- 



e 2 ,e 



'1 > 



-1 J ' 



(BIO) 



(Bll) 



(B12) 



n=l 



n,q 1 n+3,9 



* ff,_g+g ".-^ n + 3,- g A (B9) 



The unit vectors e" and e are 2d vectors so the .B vectors 
are 6 dimensional. The lattice constant of the kagome 
lattice is twice the bond length so it is 2a. 

The central force part of this dynamical matrix is 
the same as the nearest neighbor part of the dynami- 
cal matrix in Ref. 43[ as well as the central-force part in 
Ref. 0. 

The Fourier transform on the kagome lattice follows 
similar rules defined for the hexagonal lattice, albeit with 
a different unit cell area vq f- = 2\/3a 2 • 



[1] D. A. Fletcher and D. Mullins, Nature 463, 485 (2010). 
[2] S. Weiner and H. D. Wagner, Annu. Rev. Mater. Sci. 28, 
271 (1998). 

[3] M. E. Davis and R. F. Lobo, Chem. Mater. 4, 756 (1992). 
[4] P. D. Yurchenco and J. C. Schittny, FASEB J. 4, 1577 
(1990). 

[5] K. D. Hammonds, M. T. Dove, A. P. Giddy, V. Heine, 
and W. B., Am. Mineral. 81, 1057 (1996). 

[6] T. A. Witten and L. M. Sander, Phys. Rev. Lett. 47, 
1400 (1981). 

[7] D. W. Schaefer and K. D. Keefer, Phys. Rev. Lett. 56, 
2199 (1986). 

[8] S. Horike, S. Shimomura, and S. Kitagawa, Nat. Chem. 
1, 695 (2009). 



[9] M. E. Davis, Nature 417, 813 (2002). 
[10] J. Lee, O. K. Farha, J. Roberts, K. A. Scheidt, S. T. 
Nguyen, and J. T. Hupp, Chem. Soc. Rev. 38, 1450 
(2009). 

[11] T. A. Schaedler, A. J. Jacobsen, A. Torrents, A. E. 

Sorensen, J. Lian, J. R. Greer, L. Valdevit, and W. B. 

Carter, Science 334, 962 (2011). 
[12] K. Y. Lee and D. J. Mooney, Chem. Rev. 101, 1869 

(2001). 

[13] K. Sun, A. Souslov, X. Mao, and T. C. Lubensky, Proc. 

Natl. Acad. Sci. U. S. A. 109, 12369 (2012). 
[14] V. Kapko, C. Dawson, M. M. J. Treacy, and M. F. 

Thorpe, Phys. Chem. Chem. Phys. 12, 8531 (2010). 
[15] A. Sartbaeva, S. A. Wells, M. M. J. Treacy, and M. F. 



16 



Thorpe, Nat. Mater. 5, 962 (2006). [32 
G. N. Greaves, A. L. Greer, R. S. Lakes, and T. Rouxel, [33 
Nat. Mater. 10, 823 (2011). [34 
L. M. McDowellboyer, J. R. Hunt, and N. Sitar, Water 
Resour. Res. 22, 1901 (1986). [35 

G. Ernst, C. Broholm, G. R. Kowach, and A. P. Ramirez, [36 
Nature 396, 147 (1998). 

J. N. Grima, R. Gatt, V. Zammit, J. J. Williams, K. E. 
Evans, A. Alderson, and R. I. Walton, J. Appl. Phys. [37 
101 (2007). 

J. N. Grima, R. Jackson, A. Alderson, and K. E. Evans, [38 
Advanced Materials 12, 1912 (2000). [39 
J. C. Maxwell, Philos. Mag. 27, 294 (1864). [40 
J. D. Joannopoulos, P. R. Villeneuve, and S. H. Fan, 
Nature 386, 143 (1997). [41 
A. Moroz, Phys. Rev. B 66, 115109 (2002). 

A. J. Garcia- Adeva, New J. Phys. 8, 86 (2006). [42 
J. F. Galisteo-Lopez, M. Ibisate, R. Sapienza, L. S. 
Froufe-Perez, A. Blanco, and C. Lopez, Adv. Mater. 23, [43 
30 (2011). 

Q. Chen, S. Bae, and S. Granick, Nature 469, 381 (2011). [44 

H. Cohn and A. Kumar, Proc. Natl. Acad. Sci. U.S.A. [45 
106, 9570 (2009). 

E. Edlund, O. Lindgren, and M. N. Jacobi, Phys. Rev. [46 
Lett. 107, 085501 (2011). 

S. Torquato, Soft Matter 5, 1157 (2009). [47 
K. P. Velikov, C. G. Christova, R. P. A. Dullens, and 
A. van Blaaderen, Science 296, 106 (2002). [ 
[31] S. C. Glotzer and M. J. Solomon, Nat. Mater. 6, 557 
(2007). 



F. Romano and F. Sciortino, Nat. Commun. 3 (2012). 
A. V. Tkachenko, Phys. Rev. Lett. 89, 148303 (2002). 
X. Mao, Q. Chen, and S. Granick, Nature Materials 7, 
217 (2013). 

P. G. De Gennes, Rev. Mod. Phys. 64, 645 (1992). 

S. Jiang, Q. Chen, M. Tripathy, E. Luijten, K. S. 

Schweizer, and S. Granick, Advanced Materials 22, 1060 

(2010) , times Cited: 94 Si. 

J. Villain, R. Bidaux, J. -P. Carton, and R. Conte, J. 

Phys. France 41, 1263 (1980). 

C. Henley, J. Appl. Phys. 61, 3962 (1987). 

A. Chubukov, Phys. Rev. Lett. 69, 832 (1992). 

J. N. Reimers and A. J. Berlinsky, Phys. Rev. B 48, 9539 

(1993). 

A. H. Castro Neto, P. Pujol, and E. Fradkin, Phys. Rev. 
B 74, 024302 (2006). 

Q. Chen, E. Diesel, J. K. Whitmer, S. C. Bae, E. Luijten, 
and S. Granick, J. Am. Chem. Soc. 133, 7725 (2011). 
X. Mao and T. C. Lubensky, Phys. Rev. E 83, 011111 

(2011) . 

P. N. Pusey and W. Van Megen, Nature 320, 340 (1986). 
X. Mao, O. Stenull, and T. C. Lubensky, Manuscript in 
preparation (2011). 

T. C. Lubensky and P. Chaikin, Principles of condensed 

matter physics (Cambridge University Press, 2000). 

C. P. Broedersz, X. Mao, T. C. Lubensky, and F. C. 

MacKintosh, Nat. Phys. 7, 983 (2011). 

X. Mao, O. Stenull, and T. C. Lubensky, Manuscript in 

preparation (2013). 



